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»S-  ABSTRACT 


A  method  is  developed  for  remotely  determining  the  average  transverse 
wind  velocity  and  the  atmospheric  structure  constant  {strength  of  turbulence) 
at  N  points  along  a  line-of-sight  path.  The  technique  avoids  the  basic 
instability  problem  that  was  encountered  in  previous  work,  limiting  the  cal¬ 
culations  to  one  or  two  points.  Linear  integral  equations  relate  the  data, 
the  amplitude  correlation  function  and  the  amplitude  and  phase  structure 
functions,  with  the  unknown  structure  constant  and  wind  velocity.  The  stand¬ 
ard  inversion  method  leads  to  large  variations  in  the  unknown  for  small  data 
errors;  thus,  the  problem  is  ill  posed.  To  counteract  this,  a  statistical 
inversion  procedure  is  developed  that  is  dependent  upon  a  priori  knowledge  of 
the  statistics  of  the  unknowns.  The  error  in  the  final  solution  can  also  be 
predicted  by  computer  simulation.  For  example,  with  an  input  error  of  one 
percent,  the  RMS  error  in  the  unknown  will  be  chi  the  order  of  ten  percent. 
This  is  an  increase  in  accuracy  of  ten  orders  of  magnitude  over  the  standard 
inverse  moment  method. 
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I.  Introduction 

It  is  proposed  that  the  atmospheric  structure  constant  and  the  average 
transverse  wind  velocity  be  determined  remotely  through  line  of  slight  micro¬ 
wave  or  optical  scattering  measurements  [1].  The  method  will  yield  values 
of  the  unknown  at  any  number  of  points,  along  the  path,  between  the  source 
and  detector-  This  is  opposed  to  previous  methods  where  the  atmospheric 
parameters  are  predicted  for  ooly  one  or  two  points  [2,3].  The  method  has 
amplication  in  the  detection  of  clear  air  turbulence ,  the  study  of  the  basic 
atmospheric  turbulence  properties  and  in  the  everyday  measurement  of  commaro 
neteorological  parameters-  The  measurement  system  is  shown  in  Fig-  1.  The 
source  which  is  located  at  the  origin  produces  a  wave  propagation  in  the 
x-direction.  The  beam  can  take  the  form  of  a  plane,  spherical  or  beam  wave. 
The  detecting  array  is  located  at  some  distance  L  from  the  source.  It  con¬ 
sists  of  a  set  of  point  receivers  in  a  horizontal  array,  perpendicular  to 
the  x-axis.  The  signal  incident  on  the  different  array  elements  can  be 
correlated  spatially  or  temporally  to  obtain  the  statistical  properties  cf 
the  scattered  signal.  In  the  general  case,  the  statistics  are  expressed  in 
terms  of  the  structure  function,  i.e. , 

Df(?i,£2,ti.,t2)  =  E{|fCri,tj)  -  f(r2,t2)|2}  (1) 

where  r  and  t  denote  the  spatial  and  temporal  coordinates,  and  E  denotes  the 
expected  value  and  f  is  the  random  quantity.  This  function  is  useful  in  that 
it  is  expressed  roly  in  terms  of  lri-r2)  and  (t^-t2)  when  the  atmosphere  is 
assumed  to  be  locally  stationary  (locally  homogeneous ) ,  the  usual  case.  The 
more  conaro  correlation  function  expressed  as 

Bf(r1-r2,t1-t2)  =  Eiftrjjtj)  f*(r2,t2)>  (2) 

is  also  used.  It  is  a  function  of  (ri-r2)  and  (ti~t2),  and  is  realted  to  the 
structure  function  by  Eq.  (3)  if  the  atmosphere  is  assumed  to  be  strictly 
stationary  and  the  incident  beam  takes  the  form  of  a  plane  or  spherical  wave. 

Df(ri-r2,ti-t2)  =  2  B^OjO)  -  2  Bf(ri-r2,tj-t2)  .  (3) 


1 


2 


The  statistics  of  the  scattered  bean  can  be  related  to  three  basic 

atmospheric  parameters.  They  are  the  atmospheric  structure  constant,  C  2(x), 

the  atmospheric  turbulence  spectrum,  »n0/(ic) ,  and  -be  average  transverse  wind 

velocity,  V(x).  C^2  is  actually  the  square  of  the  structure  constant;  however, 

for  convenience  it  is  simply  referred  to  as  the  structure  constant.  C  2(x), 

n 

which  represents  the  strength  of  turbulence  along  the  path  and  V(x)  are  both 
considered  to  be  smoothly  varying  along  the  transmission  path;  they  are  the 
quantities  to  be  determined  by  the  inversion  method.  +^°\k)  is  assumed  to 
take  the  form  of  the  Kolmogorov  spectrum;  i.e. , 


♦<»>  .  «-‘1/3  . 


(4) 


This  simple  spectral  form  is  required  in  order  to  obtain  integral  equations 
in  closed  form.  The  turbulence  parameters  came  about  through  the  assumption 
that  the  atmospheric  index  of  refraction  has  a  slight  random  variation  about 
its  mean  value.  This  can  be  expressed  as 


n(r,t)  =  1  +  n;(r,t) 


(5) 


where  n  is  the  index  of  refraction  with  average  one  and  nj  the  random  variation. 
It  is  assumed  that 


|nij  «  1  . 


(6) 


This  and  the  "frozen-in”  hypothesis  are  used  in  the  derivation  of  the  integral 
equations  that  relate  the  atmospheric  parameters  to  the  beam  parameters. 


II.  Integral  Equations 

The  relationship  between  the  beam  parameters  and  the  atmospheric  parameters 
for  plane  and  spherical  waves  was  originally  developed  by  Tatarski  [4]  and 
extended  to  include  the  case  of  an  incident  beam  wave  by  Ishiaaru  [5].  The 
method  commences  with  the  development  of  the  wave  equation  where  the  index  of 
refraction  is  in  the  form  of  Eq,  (5).  Rytov’s  method  of  small  perturbations 


3 

is  used  to  separate  the  incident  from  the  scattered  f5eld.  The  time  dependence 
of  the  index  of  refraction  is  removed  through  the  "frozen- in"  hypothesis,  and 
the  solution  is  obtained  by  use  of  a  spectral  technique.  For  the  general  case 
of  an  incident  beam  wave,  the  correlation  functions  and  the  structure  functions 
become 


BA  }  =  V2  0-033  |  dn  Cn2(n)  |  KdK  {U0(«cF>  +  Jo(kP*)J  |h|2 
0  0 

±  J  (kQ)H  ±  J  (kQ*)h“}  *(o)(r) 
o  o  n 


(7) 


}  = 


L  • 

2i2  0.033  j  dn  C  2(n)  (  xdic  {[J  (i2y.Kpj)  +  J  (i2y.icp2)  -  J  (kP) 

J  nj  o  i  o  i  o 

0  0 

(8) 

-  jo(k?*)]  | H| 2  ±  [1  -  Jq(kQ)3  H2  ±  fl  -  jo(kq*)]  h*2]  ♦^c)(-) 


where  the  upper  sign  is  for  amplitude  data  and  the  lower  for  phase  data. 
The  variables  in  £qs.(7)  and  (8)  are 


P  = 


=  (yjd  +  iYiyc  +  V  T)2  +  (yrzd  +  iy.zc  +  v2t)2 


(9) 


Q2  =  (yyd  +  VyT)2  +  (yzd  +  VzT)2 


y*  -  yi  -  y2 


yc  =  yi  +  ya 


z,  =  Zi  -  z2 


zc  =  Z!  +  z2 


Pi  =  (yi2  +  y2 2)3s 


p2  =  (zi2  +  Z22)35 


1 H j 2  =  k2  exp(Ti  ^  k2} 


H2  =  -k  exp  (iy  ~  r2) 


i  +  ion 

Y=m5 L  =  Yr  +  iYi 


Yr  = 


1  -  02L  +  [(oj2  +  022)  L  -  02T,n 
(1  -  02L)2  +  (ojL)2 


ai(L-n) 


Yi  =  - 


(1  -  02 L)*-  +  (ajL)2 


“1  = 


vW 


°2  =  r- 


V  and  V  are  the  components  of  the  average  wind  velocity,  the  y's  and  z’s 

y  z 

are  the  receiver  coordinates,  k  is  the  atmospheric  wave  number,  k  is  the 
electromagnetic  wave  number,  Wq  is  the  radius  of  the  transmitting  aperture, 
and  Kq  is  the  position  of  the  bean  focus. 

If  $^°\»c) ,  the  atmospheric  spectrum,  is  taken  as  in  Eq.  (4),  the  following 
integral  may  be  used: 


v 

! 


kV  exp(-OK-)  J  (8k)  d»: 

o 


r( 


u  +  li 


U  +  1 


2  a 


(10) 


where 


Re(a)  >  0 


and  Re(u)  >  -1. 


Notice  that  while  the  equations  do  not  satisfy  the  conditi  ,n  Re(y)  >  -1,  the 
integrands  approach!  zero  for  small  k.  (The  phase  carrel  .ion  function  is  an 
exception,  and  will  not  converge  tinder  this  assumption. ;  Further  (tsservation 
reveals  that  the  integrand  can  be  analytically  continued  to  v  =  3/3.  After 

completing  the  integration  over  the  atmospheric  spectrum  the  ccrrelatico  and 
structure  functions  become 


5 


*  0.033  k2»2r(-5/6)L  Be  [  to  C  3(x)  IjFi  [-5/6  ;1; - — - ! 

J0 


T  5/6  02  L  5/6, 

[-6.  £(l-x)j  -jFil-S/ejl; - S - ]  U6  £  (l-x)j  }  (11) 

Xk  416  |  (1-x)  * 


^  p  2 

'(ri,r2,T)  =  n.033  k^rf-S/BjL  Re  j  dx  Cc2(x)  [{jlF^-5/6;  1;  -  — ] 

a  n  V  - 1-x) 


f  iF li-5/6  ;  1;  - 


°iP2 


p2  i  5/6 

]  -  21F1[-5/6  ;  1  ; - f- - ]}  [-6.  £  (1-x)] 

■l*  «•  -  *  i.  r  "  r.  — \  i  * 


f  (1-*> 


**i  -k-  (1-x) 


5/6-1 


?  2{l  -  iFit-5/6  ;  1  ; - f - ]}  [i6  £  (1-x)] 

4i6  £-  (1-x)  *  J 


(12) 


where  the  path  length  has  been  normalized  to  vary  from  zero  to  one,  and 


1  *  ioLx 
1  +  i  o  L 


(13) 


As  can  be  seen  from  Eq.  (11)  and  (12)  the  unknown  structure  constant  and 
wind  velocity  are  contained  within  integral  equations;  thus,  the  solution 
calls  for-  an  inversion  procedure.  Upon  closer  inspection  it  is  seen  that  this 
is  a  formidable  task.  Both  unknowns  are  involved  in  the  integral  equations, 
the  structure  constant,  C^2 ,  in  a  linear  fashion  and  the  wind  velocity,  V,  in 
a  non-linear  one.  To  simplify  the  equations  it  would  be  desirable  to  find 
them  in  terms  of  one  unknown  or  the  other,  and  both  in  a  linear  fonn.  This 
can  be  done  by  noting  that  the  wind  velocity  always  occurs  in  conjunction  with 
the  time  delay  variable,  t;  hence  the  structure  constant  can  be  found  independ¬ 
ent  of  the  wind  velocity  by  taking  t  to  be  zero.  This  results  in  a  change  in 
only  the  P  and  Q  variables  in  Eq.  (9).  They  become  Pq  and  Qp  respectively. 


»& 


F.  •  ! 


The  linearization  of  the  wind  velocity  is  accomplished  through  a  method  proposed 
by  Shen  [6].  First  differentiate  Eqs.  (7)  and  (8)  with  respect  to  t.  The  equa¬ 
tion  is  then  linearized  by  eliminating  the  V(x)  term  from  the  kernel.  This  is 


done  by  equating  t  to  zero;  the  result  is  a  linear  equation  in  Cq2(x)  V(x). 


3  A  -  . 

37=,  <r‘-l'2-T) 


=  0.033  k2*2?  (1/6)  L 


(14) 


T=0 


He 


/  Cn2(x) 


V(x) 


i  P  2 

{iFilr;2;-  0 


4*i  £  (l-x) 


3  l-&i  f  d-x)j 


-1/5 


l*rra  +  i4ircJ 


T1 


4i6  £  (l-x) 


L  - 1/6  , 

-]  (16  £  (l-x)]  (fir.)} 


where 


rd  =  <yi  -  y2>j  ♦  (zi  -  ^2)k 
rc  =  <yi  +  >’2^3  +  (21  +  z2)k 

Po2  =  (Vd  +  iTiyc)Z  +  (Vd  +  iYiZc)2 
Q02  =  (yyd)2  +  (yyc)2 


(15) 


While  cn2(x)  and  V(x)  cannot  be  obtained  independently  in  this  case,  the 


wind  velocity  can  be  determined  since  C^2  is  already  known  from  the  previous 


equations.  In  this  derivative  format,  it  is  found  that  the  correlation  function 
and  the  structure  function  are  related  in  a  very  simple  way,  i.e.. 


37  Vri’r2,T) 


*  -  2  37  Vri,r2,T) 

T=0  T=0 


(16) 


Thus,  either  one  can  be  used  in  the  solution  for  the  wind  velocity. 


III.  The  Inversion  Method 

The  equations  describing  the  scattering  of  waves  in  the  atmosphere,  as 
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derived  in  the  last  section,  are  found  to  take  the  form  of  a  Fredholm  integral 
equation  of  the  first  kind.  The  general  form  of  this  integral  equation  is 

B 

g(y)  =  j  K(x,y)  f(x)  dx,  o’  f  y  f  6'.  (17) 

a 

g(y)  is  a  known  function  or  data,  K(x,y)  is  the  kernel  of  the  integral  equation, 
f(x)  is  the  unknown  and  a,  B,  a* ,  B*  are  fixed  constants.  The  equation  ccm  be 
solved  analytically  if  the  kernel  is  very  simple  or  if  it  can  be  expressed  as  a 
complete  set  of  orthogonal  functions.  When  the  kernel  is  more  complicated  the 
use  of  numerical  methods  is  usually  necessary,  and  the  moment  method  is  commonly 
employed.  This  method  is  developed  by  expanding  the  integral  equation  into  N 
simultaneous  equations  in  N  unknowns,  and  contracting  into  matrix  form,  as  shown. 


OQ 

II 

(18) 

where 

g  =  tg(yi)l 

(19) 

f  =  [f(x.)J 

and 

A  =  [W(x.)  K(xiSy.)] 

where  the  braces  enclose  tr.e  elements  of  the  g,  f,  and  A  matrices.  W(x^)  is  a 
weighting  fui'ct’on  dependent  upon  the  quadrature  expansion  of  the  integral. 

In  subsequent  steps  W(x. )  will  not  be  shown  since  it  can  be  carried  with  the 
kernel.  ’:'he  x^’s  and  y^'s  are  discrete  values  in  the  range 

a  5  x.  f  B  (20) 

and 

a1  5  y.  »  B’ . 

i 

The  solution  of  Eq.  (18)  can  be  easixy  obtained  by  numerically  inverting  the  A 
matrix;  however,  it  is  soon  discovered  that  the  results  are  highly  unstable, 
and  do  not  represent  the  unknown  by  any  stretch  of  the  imagination.  The  problems 
arise  from  the  errors  associated  with  the  data  and  those  introduced  in  the  quad¬ 
rature  expansion  and  the  inversion  process.  This  can  be  described  mathematically 


8 


by  the  following  uniqueness  argument.  Suppose  that  two  sets  of  data,  gj  and  g2, 
correspond  to  two  sets  of  unknowns ,  fj  and  f 2,  where  f2  =  f 1  +  M  sin  («x).  In 
integral  equation  form 


and 


8 

gl(y)  =  |  K(x,y)  fxCx)  dx 


8 

g2(y>  =  |  K(x,y)  f2(x)  dx  = 


8 

r 

K(x,y)  tfi(x)  +  H  sin(ttx)]  dx  . 
a 


(21) 


(22) 


Equation  (22)  can  be  expanded  to 

6  8 


g2(y)  =  |  K(x,y)  fx(x)  dx  +  |  K(x,y)  H  sin(ux)  dx 


or 


p 

g2(y)  =  gi(y)  +  |  K(x,y)  W  sin(ux)  dx  . 


(23) 


(24) 


For  any  constant  u  can  be  chosen  large  enough  so  that  the  integral  of  the 
kernel  and  the  rapidly  varying  sine  term  average  to  an  arbitrarily  small  constant, 
or 

g2(y)  =  gl(y)  +  e  .  (25) 


Thus,  for  two  sets  of  data  that  vary  by  only  some  small  experimental  error,  the 
values  of  the  unknowns  can  differ  by  W  sin(ux),  a  highly  oscillatory  function 
of  great  magnitude.  This  indicates  that  the  solutions  of  the  Fredholm  integral 
equation  of  the  first  kind  are  not  unique  when  experimental  errors  are  taken 
into  account . 

To  compensate  for  this  problem  one  must  account  for  the  errors  in  the  data 
and  the  unknown.  This  can  be  done  by  modifying  Eq.  (18),  i.e., 

g  +  e  =  A(f  +  *>.  (26) 


c  can  be  considered  the  experimental  error  in  the  data  and  (  the  resultant  errors 
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in  the  unknown.  It  is  sometimes  convenient  to  denote  g  +  e  and  f  +  5  as 


g  =  g  +  e 

(27) 

f  =  f  +  e  . 

(28) 

g  is  the  data  actually  used  in  the  determination  of  the  unknown  since  the  true 
data  and  its  error  are  inseparable.  By  the  same  token,  f  represents  the  solu¬ 
tion  that  is  obtained  from  the  inversion  process.  Of  course  it  is  desired  that 
f  approach  f  as  closely  as  possible  and  this  condition  is  attained  by  minimizing 
C- 

The  minimization  of  5  is  accomplished  by  noting  that  the  A  matrix  in  Eq. (26) 
is  a  linear  operator;  hence  its  inverse  m>'_t  also  be  linear,  and 


f  +  £  =  B(g  +  e)  (29) 

or 

(  *  -f  ♦  B(g  t  e).  (30) 


The  B  matrix  is  an 
minimization  of  £. 
is  obtained  of  the 


unknown,  linear  operator  that  is  to  be  determined  by  the 
Multiplying  the  vector  £  bv  its  transpose,  a  square  matrix 
form 


Ci  Cl  Ci  C2 

C2  Ci  C2  C2 


*1  Cn  \ 
“■  \ 


Cn  Cl 


Cn  C2 


Cjj  Cj{ 


(31) 


Through  the  minimization  of  the  diagonal  terms,  the  B  can  be  found,  the  solution 
amounts  to  a  minimum  squared  error  method.  This  is  analogous  to  results  obtained 
by  Franklin  and  others  [7,8].  (Previously,  a  deterministic  approach  had  been  taken 
by  many  people  {9,10,11,12,13].)  This  specific  form  of  the  solution  is  found 
convenient  in  the  later  determination  of  the  wind  velocity.  Expanding  Eq.  (30) 
as  indicated  above,  an  equation  in  terms  of  square  matrices  is  obtained. 


Hi*'*1 


(32) 


1C 


«T  =  ffT  -  f(Bg)T  -  (Eg)fT  -  f(Be)T  -  (Be)fT 
+  (Bg)(Be)T  +  (Be)(Bg)T  +  (Bg)(Bg)T  +  (Be)(Be)T 


with  diagonal  tezms 


hh  -  fifi '  2fi  I  Bik*k  ‘  2fi  I  BikEk 

+  2  |  Bikgk  |  Bilel  +  |  Bik*k  |  Bit£l  *  £  BikEk  |  BiiEt  • 


(33) 


To  minimize  this  equation,  it  is  differentiated  with  inspect  to  each  element 
of  the  B  matrix,  B  ,  obtaining 


jr2  *  0  =  0  -  2fi  I  4kn%  -  2fi  I  5knEk 
in  k  k 

*  2  j[  ‘kA  |  BitE»  +  2  1  BikEk  |  5tnEt 

+  2  I  4kngk  l  BilH  *  2  l  5knEk  \  BitV 


5,  is  a  matrix  with  all  terms  zero  except  the  k,  nth  which  is  one.  Notice 
kn  Hi$i 

that  all  equations  developed  from  -sr - with  i  i  m  are  zero.  Since  both  i  and 

Hi 

n  range  from  one  to  N,  a  set  of  H  *  H  simultaneous  equations  are  obtained  for 
the  solution  of  the  B  matrix.  Doing  the  summations  over  the  delta  functions 
Eq.  (34)  becomes 


0  *  -fA  -  fiEn  *  H  l  BitEt  *  E»  l  Bik«k  *  «n  l  BilH  +  S.  |  BitEl  <35) 

l  K  l  L 


which  can  be  contracted  back  into  matrix  form  to  obtain 


0  =  -fgT  -  fe*  +  BgeT  +  BegT  +  BggT  +  BeeT.  (36) 

This  equation  relates  the  B  matrix  to  the  unknown,  the  data,  and  the  error  in 
the  data.  Since  it  is  desired  to  find  B  in  terms  of  the  statistics  of  these 
quantities  it  is  necessary  to  take  an  expected  value  to  obtain  equations  in 
terms  of  the  covariances  matricies. 
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0  =  -R_  -  Rg  +  B(R  +  R  +  R  +  R  ) 
fg  fe  ge  eg  gg  ee 

B.6  -  E(«BT> 


(37) 

(38) 


To  further  simplify  Eq.  (37)  the  propagation  theorem  is  used  (i.e. ,  if  o  =  AB 
then  R  =  AR00AT  and  R0  =  R.0A^)  obtaining 

OUX  pp  ptt  pp 

0  =  ~RffAT  -  Rf£.  +  B(ARfc  +  RfgAT  +  ARffAT  +  R  )  .  (39) 

Finally  one  solves  for  B  and  subsitutes  back  into  Eq.  (29)  with  the  result: 


f  =  (RffAT  +  Rf£>  (ARffAT  +  ARfe  +  Rf£.AT  +  S^)'1  £  .  (40) 

This  is  identical  to  the  form  presented  by  Franklin.  In  practical  situations 
the  unknown  is  independent  of  the  data  error,  requiring 


and 


(41) 


f  =  RffAT  (ARffAT  +  RMr»  £ 


(42) 


which  is  the  usual  form  of  the  inversion  equation.  As  a  special  case  assume 
that  the  errors  in  the  system  are  zero.  This  implies  R££.  =  0,  f  f ,  and  £  -*■  g; 
resulting  in 


f  =  A~l  g  ,  (43) 

the  original  integral  equation  when  no  errors  are  present. 

The  predicted  accuracy  of  the  inversion  method  can  be  found  by  taking  the 
expected  value  of  Eq.  (32).  By  using  the  matrix  propagation  theorem  and  com¬ 
bining  terms,  the  following  is  obiained. 

R«  =  Rff  *  (Rf/  *  V  (ARff*T  +  Rec  +  Rtf*T  +  “fe*'1  (ARff  *  Rcf>  <«> 


R^  is  the  covariance  matrix  of  the  error  term;  the  other  terms  have  been  defined 
previously.  The  RMS  error  of  the  predicted  value  of  the  unknown  is  then 
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H  =  tN  VR££)jlS  =  *N  .1^ i2]5*  • 


Another  method  for  error  evaluation  has  been  provided  by  Franklin.  A  is 


defined  as 


lull 

| j  f  +  £  1 J  normalized  error  in  the  unknown 
j  |  e 1 1  normalized  error  in  the  data 

Ms  +  ell 


where  ] ]h| j  is  the  norm  of  h.  If  delta  is  on  the  order  of  one  then  the  error 
generated  in  the  unknown,  by  the  inversion  method,  is  about  the  same  as  the 
error  in  the  data.  Clearly,  if  delta  is  large,  then  the  unknown  error  is  much 
greater  than  the  data  error  and  the  inversion  method  is  not  satisfactory.  Since 
the  norm  is  simply  a  number,  Eq.  (46)  can  he  modified  to 

A  =  JilLL  lif..*  -.11  .  (47) 

lleji  ||f*Ci| 

This  relation  can  be  simplified  with  the  equations 

g  +  c  =  A(f  +  O  (48) 

and 

f  +  £  =  B(g  +  e).  (49) 

If  it  is  assumed  that  the  error  in  the  unknown  is  strictly  due  to  the  errors  in 
the  data,  then  the  second  equation  will  reduce  to 

£  =  Be  (50) 

for  no  input.  In  this  case  Eq.(47)  becomes 


!*<*♦  *;i 

Ilf  *  ill 


For  the  worst  case 
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-  »«  1IWI1  (ll»<f  »  on 

-  max  i  ,  '  , * i  ®ax  \  « s  _  . ,  J 


■11*11' 


Ilf  ♦  Cll 


By  use  of  the  Schwartz  inequality,  i.e.. 


1M1 

l|a|| 


l°l  1  llAii  _ 

"Tull "  '  l,A!l  ' 


TUTr* 


Eq.(51)  can  be  simplified  further.  The  final  form  is 

AaaX  =  I  lAl  i  I lBll-  <54) 

By  defining  Ae  and  A£  to  be  the  average  error  in  the  data  and  unknown  respectively 
and  using  the  definition  of  A,  Eq. (46),  it  can  be  seen  that 

At  =  I  Ml  i|B||  Ae  .  (55) 

These  equations  lead  to  the  solutions  of  the  atmospheric  structure  constant  at 
several  points  along  the  path.  An  estimate  of  the  error  in  the  solution  is  also 
obtained. 

When  the  solution  of  the  wind  velocity  is  desired  Eq.  (14)  must  be  solved. 

Unlike  the  previous  case,  the  unknown  wind  velocity  is  associated  with  another 

variable,  the  structure  constant.  Consequently,  the  solution  of  the  equation 

will  not  yield  the  wind  velocity  profiles  directly.  Since  C^~  is  known  from  the 

previous  developments  it  should  be  possible  to  evaluate  the  wind  velocity  itself. 

It  was  found  that  the  most  accurate  method  of  solution  was  to  associate  C  2  with 

n 

the  kernel  matrix  obtaining 

g*  =  A* V  (56) 

where  g*  is  the  derivative  data,  V  is  the  unknown  and  A’  contains  the  kernel, 
the  weighting  function  and  the  structure  constant,  i.e.. 


A’  =  (W(x.)  Cn2(yj)  X‘(pi,xj)J. 


In  the  implementation  of  the  solution  it  is  again  necessary  to  acount  for  the 
errors  in  the  data  and  the  unknown;  thus. 
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g*  +  e'  =  A' (V  +  c),  (58) 

where  e'  is  the  error  in  the  derivatives  and  c  the  errors  in  the  wind  velocity. 

It  should  also  be  remembered  that  the  value  of  C  2  used  in  the  kernel  is  not 

n 

exact  but  contains  sente  error;  this  implies  that  the  A  matrix  is  actually  of 
the  form 

A’  =  [(C  2j  +  i .)  Ki.J  .  (59) 

11  J  ^  J 

(W(Xj)  is  contained  in  the  K*  term.) 

While  this  presents  more  complications  it  is  assumed  that  the  solution  is  in 
terms  of  a  linear  operator,  thus 

V  +  ?  =  B'(g’  +  e')  .  (60) 

The  coefficient  matrix,  B' ,  is  found  by  minimizing  £,  thereby  obtaining 

Vg'T  =  B’g»g,T  +  Be*E,T.  (61) 

The  procedure  is  identical  to  the  last  section  where  the  errors  in  g  are  assumed 
to  be  independent  of  V.  With  the  proper  substitutions  g  can  be  eliminated 
obtaining  an  equation  of  the  form 

WTA,T  =  B'A*WTA,T  +  B’e'e,T  (62) 

with  substitution  of  Eq. (59)  one  obtains  the  following: 

WT[(C2  +  £.)  K!  -  jT  =  B»l(C2  +  g  )  K!  ]  WT[(C2  +  g.)  K!.]T  +  B'e'e'T.  (63) 
**J  J  XJ  nJ  j  ij  ’j  a] 

The  expected  value  can  now  be  taken  with  the  reasonable  assumption  that  e,  the 
error  in  the  structure  constant  is  independent  of  the  wind  velocity. 

WClj  Ki/  +  Kwl«y*i/  -  K!.)  iytCj.  K!.]1 

+  B'Icnj  Kij3  Kij]T  +  Kij3  V^n  Kij]T  (64) 

+  B’Re,e,  +  B*E[(£5K!.)  WT(g.W.)T] 


^-SS*S£s53^« 
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It  can  be  shown  that  E(£^ )  =  0  if  the  expected  value  of  the  error  .In  the 


structure  function  is  zero.  It  is  also  assumed  that  the  errors  in  the  structure 
constant  are  related  only  to  the  errors  in  the  structure  function.  With  this 
and  with  simplification  of  the  last  term  in  Eq.  (64)  the  equation  becomes 


WCnjKij,T  =  B'lCnjKlj)  V'^i/  *  B’(V  '  «'rwk'T>  *  B’\'E- 


nj  13 


U3  13 


(65) 


where  the  dot  denotes  the  matrix  operation  defined  by 


a  •  6  =  (a..  B-.) 

13  13 


i»3  =  1*2,  •••,  N. 


(66) 


Solving  forB'  and  'ubstituting  into  Eq.  (60)  V  can  be  found 


5  ■  RVV<CnjKij)T  £<C»iKij>  ’W<5jKij>1 


+  (Btt)  *  (K'R^K*1)  +  Re,etJ-1  DCPi,x) 


(67) 


T-0 


This  expression  is  similar  to  that  derived  by  Franklin’s  method.  The  present 

inversion  method  introduces  an  additional  term  to  account  for  the  errors  in  C  2« 

n 

The  preceding  solution  gives  the  best  mean  squared  estimate  of  the  unknown 

structure  constant  and  average  transverse  wind  velocity.  The  solutions  are 

dependent  upon  the  input  data,  the  correlation  and  structure  functions  of  the 

scattered  beam  that  propagates  through,  the  medium  to  be  remotely  sensed.  The 

solutions  are  also  dependent  upon  the  statistics  of  the  unknowns  and  the  error. 

These  are  represented  by  the  covariance  matricies  R-,.  and  R  .  It  is  next  nec- 

rr  ee 

essaiy  to  determine  the  form  of  these  matricies. 


IV.  Statistical  Quantities 


In  the  development  of  the  unknown  covariance  matricies  it  is  convenient 
to  represent  the  unknown  functions  in  terms  of  a  random  Fourier  series. 


N 


f(x)  =  f  +  J  bn^an  cos(nitx)  +  a'  sin(nxx)] 
n=o 


(68) 
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where  f  is  the  mean  value  of  f(x),  the  b  's  are  fixed  constants  relating  the 

n 

magnitude  of  the  variations  to  the  mean  value,  a^  and  an’  are  independent 
random  coefficients  with  zero  mean  and  variance  one  giving  the  necessary 
variability  to  the  unknown,  and  N  limits  the  rate  at  which  the  fluctuations 
occur.  The  function  also  has  the  advantage  of  being  easily  generated  in  com¬ 
puter  simulation  schemes. 

To  find  the  covariance  matrix  one  employes  Eq.  (38)  obtaining 


where 


and 


Eff  * f2 


N 

l 


n=o 


b  2  C06  n*(x„  -  x.) 
n  13 


E(a  a  )  = 
n  m 


mn 


E(a  )  =  0 
n 


(69) 

(70) 

(71) 


as  stated  previously.  The  f  and  b^  coefficients  are  chosen  from  experience  and 
certain  realizability  conditions,  the  specifics  of  which  be  given  in  the  next 
section. 

The  covariance  matricies  of  the  error  terms  are  developed  with  the  assump¬ 
tion  that  the  e^’s  are  independent  of  each  other  yielding 


R 

ee 


<*e2  V' 


(72) 


The  covariance  matrix  of  the  error  in  the  derivative  function  is  developed  in 
a  slightly  more  general  way.  Since 


!i  =  li*  Ei 


the  derivative  is  represented  by 


gi+k  "  gi-k  .  ci+k  “  ei-k 

- +  -  ; 

X.  ,  -  X.  ,  X.  ,  -  X-  . 

i+k  i-k  l+k  i-x 


(73) 


(74) 


thus ,  e!  corresponds  to  the  last  term  of  Eq.  (74);  i.e.. 


C  •  »  *  C  •  . 

,  _  i+k _ i-k 


e:  = 

1  x. -  x. 


i+k 


i-k 


(75) 


R  ,  ,  =  -  v-  (26. .  -  5 .  .  ,  ,  6.  .  .  ,  ) 

e'e '  4R2h2  13  i+k,x-k:  i-k,.\+k 


(76) 


wr^re  h  represents  the  increment  between  the  x^'s ,  the  i  and  j  subscripts 
represent  the  i^jth  matrix  element.  The  k  represents  the  increment  over  which 
the  derivatives  of  the  data  are  taken.  If  k  is  taken  very  small  the  errors  in 
the  structure  functions  will  produce  very  large  errors  in  the  derivative  data. 
If  k  is  too  large  the  numerical  derivatives  will  not  represent  the  slope  of  the 
function,  k  is  finally  chosen  as  a  compromise  value  that  best  fits  the  partic¬ 
ular  numerical  solution. 


V.  Numerical  Evaluation 

The  numerical  evaluation  of  the  theory  is  accomplished  through  both  a  com¬ 
puter  simulation  and  the  evaluation  of  experimental  results .  For  convenience, 
the  parameters  chosen  are  modeled  after  an  operational  system  under  the  direc¬ 
tion  of  A.  T.  Waterman  at  Stanford  University.  The  system  consists  of  a  trans¬ 
mitter  located  on  the  east  side  of  San  Francisco  Bay  at  an  elevation  of  300 
meters.  To  implement  the  theory,  the  transmitter  is  assumed  to  be  located  at 
the  origin  of  a  coordinate  system  with  the  beam  propagating  along  the  x-axis. 

The  beam  shape  approximates  that  of  a  spherical  wave.  A  cross  section  of  the 
transmission  path  is  shown  in  Fig.  2.  It  is  28  km  in  length,  traversing  San 
Francisco  Bay,  with  the  terminus  era  the  west  side  of  the  bay  at  an  elevation 
of  120  meters.  The  transmission  path  is  perpendicular  to  the  longitudinal  axis 
of  the  bay,  thus  it  is  likely  that  any  wind  will  blow  transverse  to  the  path 
of  the  beam.  The  wavelength  of  the  incident  beam  is  8.6  mm;  the  transmitting 
antenna  is  a  1.5  m  diameter  paraboloid  with  a  0.4  degree  beam.  Using  these 
conditions  one  can  find  e  in  the  beam  wave  equation.  It  is  found  to  be  4.86  x 
10~3.  This  determines  the  point  at  which  t  e  incident  wave  makes  the  transition 
from  a  plane  wave  to  a  spherical  wave;  it  is  about  206  meters  from  the  transmitter. 
Since  the  path  is  28  km  in  length,  it  is  seen  that  the  incident  field  behaves  as 
a  spherical  wave  over  virtually  the  entire  path.  The  receiving  system  consists 
of  an  array  of  eight  point  receivers  located  near  the  axis  of  the  beam.  The 
array  is  perpendicular  to  the  path  and  positioned  horizontally  above  the  ground. 


>coag 
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with  element  spacings  of  400  wavelengths.  See  Fig.  3.  If  the  first  element 
is  assumed  to  be  on  the  x-axis ,  the  kf.mel  of  Eq.  (11)  takes  the  form 


K(x,p)  =  w2  0.033  k2  T(-5/6)  L 

r  p2  ,  .  r  ,5/6 

Re  xFii-5'6  ;  1; - £ - }  J-5.  £  (l-x)} 

46.  £  (l-x)  *  1  K 


(77) 


2 

-  iFii-5/6  ;  1  j - =2 - }  fi«  £(l-x)} 

4*6  ~  (l-x)  K 


5/6 


wnere 


P2  =  02  =  6  0.086  •  400  (j-1)  j  =  1,2,  •••  8 
L  =  28  x  io2 
o  =  4.86  x  10”3 


(78) 


k  =  2W/0.086 

j  corresponds  to  the  position  of  the  jth  element  in  the  receiving  array.  The 
kernel  was  evaluated  numerically  to  an  accuracy  of  8  decimal  places.  Both 
ascending  and  asymptotic  series  are  used  to  evaluate  the  hypergeometric  func¬ 
tion,  the  choice  depending  upon  the  magnitude  of  the  argument. 

After  evaluation  of  the  kernel,  the  atmospheric  structure  constant,  CQ2, 
can  be  determined  from  the  matrix  equation 


[BA(p.)l  =  [A(p.,x.)l  (Cn2(x.)]  (79) 

through  use  of  the  statistical  inversion  method.  The  A  matrix  is  the  combina¬ 
tion  of  the  kernel  and  the  weighting  function.  It  wi  11  be  shown  that  when  this 
inversion  procedure  is  used,  the  predicted  errors  in  the  unknown  are  about  ten 
times  the  error  in  the  data.  This  is  referred  to  as  -the  sensitivity  of  the 
inversion  method.  For  comparison,  it  would  be  interesting  to  know  the  sensitivity 

if  a  standard  matrix  inversion  were  used  to  determine  C  2.  This  can  be  found 

n 

through  the  product  of  the  matrix  norms,  as  in  Eq.  (55).  In  this  case,  the  B 
matrix  is  sinply  A”1,  since  in  the  standard  matrix  inversion,  it  is  assumed  that 
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f  =  A-1  g.  (80) 

The  result  of  the  computation  is  ) | A{ }  | | A" * | j  =  7.44  x  10 11  indicating  that 

for  a  one  percent  data  error,  the  error  in  the  unknown  will  be  on  the  order  of 
10 11  percent.  This  certainly  leaves  some  doubt  as  to  the  existence  of  any  method 
that  could  counteract  such  large  instabilities.  Consequently,  to  confirm  the 
usefulness  of  the  inversion  method,  a  computer  simulation  was  implemented. 

To  reasonably  evaluate  the  inversion  method,  the  true  value  of  the  unknown 
structure  constant  must  be  known  at  many  points.  This  true  value  was  generated 
through  the  use  of  an  equation  similar  to  Eq.  (68);  i.e., 

1 

C  2(x)  =  10"14  +  0.4  x  io* 14  J  [a  cos(nirx)  +  a^  sin(nnx)}  (81) 

n  n=o 

where  x,  the  normalized  path  length,  varies  from  zero  to  one.  The  a  and  a' 

5i  n 

are  computer  generated,  Gaussian  random  numbers  with  zero  mean  and  variance  one. 
The  10“ 14  represents  a  typical  average  value  of  Cn2  as  derived  from  atmospheric 
experiments.  After  C^2  is  generated  at  eighteen  equispaced  points  along  the 
path,  the  true  data,  3A(p),  is  calculated  by  matrix  methods  from  Eq.  (79).  To 
model  the  errors  that  are  inherently  present  in  any  real  measurement  device,  the 
data  are  perturbed  by  adding  a  certain  amount  of  error  to  them.  The  magnitude 
of  the  error  is  defined  as  the  standard  deviation  of  the  error  deviated  by  the 
peak  of  the  correlation  curve.  It  should  be  evident  that  the  errors  between  0.1 
and  1.0  percent  would  be  typical  in  experimental  cases;  this  leads  to  errors  on 
the  order  of  tens  of  percents  in  the  tail  of  the  curve.  The  perturbed  data, 
which  approximates  the  true  data,  are  used  in  the  inversion  method  to  obtain  a 
prediction  for  Cr2;  this  prediction  is  denoted  by  C^2.  The  accuracy  of  the  pre¬ 
diction  is  found  by  comparing  C^2  with  C^2 .  The  computer  simulation  scheme  is 
diagramed  in  Fig.  4.  An  example  of  the  results  of  the  simulation,  for  the  a  jios- 
pheric  structure  constant  is  shewn  in  Figs.  5  and  6.  The  errors  involved  in  each 
case  are  indicated  in  the  figures.  For  data  errors  between  0.1  and  1.0  percent 
the  structure  constant  is  predicted  quite  accurately.  Even  for  large  errors  in 
the  correxition  function,  the  points  denoted  by  the  squares ,  the  predicted  value 
of  the  structure  constant  corresponds  to  the  mean  value  of  the  true  curve. 
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The  sensitivity  of  the  inversion  can  be  represented  in  another  way  as  shown 
in  Fig.  7.  The  average  percent  error  in  the  unknown  is  plotted  versus  the  aver¬ 
age  percent  error  in  the  data.  Three  curves  are  shown:  a)  the  error  as  pre¬ 
dicted  by  the  product  of  the  matrix  norms,  b)  the  errors  as  predicted  by  the 

covariance  matrix  of  R  ,  and  c)  the  actual  errors  as  derived  from  the  computer 

nn 

simulation.  When  the  data  error  is  between  0.1  and  1.0  percent,  the  unknown 
error  lies  between  1.0  and  10.0  for  the  worst  case  (highest  curve).  This  is  a 
remarkable  improvement  over  the  errors  generated  by  the  standard  matrix  inver¬ 
sion  method,  as  indicated  earlier  in  this  section.  An  interesting  feature  occurs 
for  the  larger  data  errors;  the  unknown  errors  seem  to  be  limiting.  This  indi¬ 
cates  that  the  even  for  large  data  errors,  the  statistical  inversion  method  is 
useful  for  predicting  the  average  value  of  the  unknown. 

The  structure  constant  can  be  found  from  data  other  than  the  correlation 
function;  namely,  the  amplitude  and  phase  structure  functions.  In  tfteory,  the 
phase  correlation  function  could  also  be  used;  however,  due  to  the  foregoing 
form  of  the  Kolmogorov  spectrum,  it  does  not  converge.  An  error  analysis  has 
been  done  for  data  in  the  form  of  the  amplitude  and  phase  structure  functions. 
Figure  8  shows  the  results,  comparing  them  with  the  correlation  function  case. 

It  is  found  r-'nat  the  unknown  errors  predicted  in  the  case  of  the  amplitude 
structure  function  are  slightly  higher  than  those  from  the  correlation  function 
expansion.  On  the  other  hand,  the  error  produced  by  the  phase  data  is  much  too 
large  to  be  of  any  use. 

The  other  parameter  to  be  measured  through  the  inversion  process  is  the 
average  transverse  wind  velocity.  To  evaluate  the  accuracy  of  the  inversion, 
an  error  analysis  similar  to  that  of  the  last  section  was  studied.  Since  the 
three  methods  of  error  prediction  (matrix  norm,  covariance  matrix,  and  computer 
simulation)  have  been  shewn  to  agree  quite  well,  only  one  will  be  used  to  eval¬ 
uate  the  wind  velocity  equations.  The  simplest  of  the  three  is  the  procedure 
employing  the  product  of  the  matrix  norms.  The  norm  of  the  A  matrix  is  obtained 
from  the  derivative  form  of  the  kernel  matrix;  the  norm  of  the  B  matrix  is  cal¬ 
culated  from  Eq.  (67).  Since  this  equation  is  rather  complicated,  it  is  sim¬ 
plified  by  assuming  C^2  constant  with  zero  error.  The  results  of  the  error 
analysis  are  shown  in  Fig.  9.  Curves  are  shewn  for  both  amplitude  and  phase 
data.  The  errors  predicted  for  the  amplitude  data  are  quite  low;  those  pre- 
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dieted  in  the  case  of  the  phase  are  quite  high.  This  reinforces  the  conclusion 
that  the  phase  data  are  not  suitable  for  the  inversion  method. 

Data  obtained  from  an  experiment  conducted  by  J.  C.  Harp  [14]  at  the 
Stanford  Electronics  Laboratories  will  now  be  evaluated  to  predict  the  wind 
velocity  and  structure  constant,  at  several  points  between  the  transmitter  and 
the  receiver.  The  data  consists  of  a  set  of  correlation  curves.  To  adapt  them 
to  the  problem  one  observes  the  value  of  t -e  correlation  function  and  its  deriv¬ 
ative  at  the  point  t=0.  In  Harp's  paper  -lata  points  are  found  for  seven  receiver 
separations.  These  data  are  interpolated  and  shown  in  Figs,  j.0  and  11.  Notice 
that  the  correlation  curves  are  normalized  to  one.  This  being  the  case,  it  is 
inpossible  to  determine  the  average  value  of  the  structure  constant  in  the  two 
cases  to  be  studied.  This  limits  one  to  examining  the  shape  and  relative  magni¬ 
tudes  of  the  struct we  constant  along  the  path.  The  average  value  of  the  wind 
velocity  is  not  lost  through  the  normalization  process.  Tc  denarmalize  the 
correlation  curves,  the  peak  of  each  is  assumed  to  have  a  typical  value  of  1.0 
x  10-2.  The  other  parameters  necessary  for  the  inversion  are  the  statistics  of 
the  unknowns  and  the  assumed  values  for  the  average  data  error.  The  covariance 
matricies  modeling  the  structure  constant  and  the  wind  velocity  are  shown  below. 

10 

Scc  =  1(T28  {l.O  +  (0.04)2  l  cos [ns(x.-x.)]}  (82) 

n=o  1  ' 

6 

Sw  =  25{l.O  +  (O.i)2  l  coslnvU.-Xj)]}  (83) 

The  average  value  of  Cn2  is  assumed  to  be  10'14;  the  (0.04)2  in  Ec.  (82)  is 
developed  from  the  constraint  that  C^2  is  greater  than  or  equal  to  zero.  The 
average  value  of  the  wind  velocity  was  taken  to  be  zero,  since  its  direction 
can  vary;  its  standard  deviation  was  taken  as  5  meters/sec.  The  magnitude  of 
the  error  in  the  correlation  function  was  taken  to  be  5  percent;  this  is  re r- en¬ 
able  considering  the  amount  of  interpolation  that  is  necessary. 

Using  the  above  parameters,  the  atmospheric  structure  constant  was  evaluated 
using  Eq.  (42).  The  results  are  shown  in  Figs.  (12a)  and  (13a).  Equation  (67) 
was  used  to  determine  the  average  transverse  wind  velocity.  These  results  are 
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shown  in  Figs.  (12b)  and  (13b).  The  structure  constant  curves  can  be  inter¬ 
preted  physically  by  observing  the  topography  over  which  toe  turbulence  was 
formed.  A  cress  section  of  the  transmission  system  was  shown  in  Fig.  (2). 

It  represents  San  Francisco  3 ay  bordered  on  both  sides  by  hills.  If  a  wind 
were  blowing  up  or  down  the  bay,  a  velocity  gradient  would  be  formed  from  the 
difference  in  the  velocity  in  the  center  of  the  bay  and  the  smaller  velocity 
that  would  occur  near  the  hills.  This  gives  rise  to  higher  turbulence  near 
the  sides  of  the  valley  than  in  the  center.  At  the  same  time,  the  turbulence 
very  near  the  sides  would  be  reduced  because  the  larger  eddies  could  not  exist 
at  that  location.  These  facts  are  reflected  in  the  curves  representing  the 
structure  constants.  The  curves  also  compare  reasonably  well  with  the  results 
developed  by  Harp.  The  structure  constant  as  determined  by  Harp  was  found  at 
three  points  along  the  path:  near  the  transmitter,  in  the  center,  and  near  the 
receiver.  The  curves  were  assumed  to  be  constant  in  these  regions.  As  can  be 
seen,  the  minimums  and  maxiEuns  of  the  structure  constant  curves ,  as  derived 
from  the  inversion  method  agree  with  those  obtained  by  Harp.  One  should  remem¬ 
ber  that,  due  to  the  normalization  of  the  correlation  curves,  the  plots  cannot 
be  compared  in  absolute  magnitude. 

The  wind  velocity,  as  predicted  by  the  statistical  inversion  method,  can 
be  closely  compared  with  that  found  by  Harp.  The  velocity  curves  are  not 
affected  by  the  normalization  of  the  correlation  function,  and  the  curves  pre¬ 
sented  reflect  both  the  general  shape  and  the  absolute  magnitude  of  the  wind 
velocity.  Before  a  comparison  of  the  curves,  it  should  be  ccntiooed  that  Harp’s 
wind  velocity  data  were  inferred  by  predicting  the  velocity  at  only  two  points 
near  the  ends  of  the  path,  and  assuming  a  smooth  variation  between  the  two 
values.  Thus,  the  velocity  plots  will  not  necessarily  compare  near  the  center 
of  the  path.  As  seen  by  inspection  of  Figs.  (12b)  and  (13b),  the  velocity  plots 
are  in  a  very  close  agreement.  Variations  do  occur  near  the  center  of  the  propa- 
g-  ion  path,  and  in  some  cases  at  other  points;  however.  Harp’s  curves  seem  to 
indicate  the  average  value  that  would  be  obtained  from  the  more  general  statis¬ 
tical  inversion  curves. 
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VI.  Conclusion 

The  value  of  the  statistical  inversion  Tsethod  for  predicting  the  atmos¬ 
pheric  structure  constant  and  the  average  transverse  wind  velocity  at  several 
points  along  the  path  has  been  demonstrated  by  computer  simulation  and  by 
application  to  data  takei  under  normal  atmospheric  conditions.  In  both  cases 
it  was  found  that  the  predicted  value  of  the  unknown  mas  within  ten  percent 
of  the  true  value  for  reasonable  data  errors.  It  was  also  determined  that  for 
larger  errors  the  predicted  solutions  correspond  to  the  average  value  of  the 
true  curve.  In  addition,  the  error  analysis  has  shown  that  the  use  of  ampli¬ 
tude  data  in  the  inversion  method  leads  to  solutions  that  are  ten  times  more 
accurate  than  Those  obtained  from  phase  fluctuation  data.  On  this  basis,  phase 
data  are  deemed  inappropriate  for  use  in  the  inversion  method.  From  these 
results  it  can  be  concluded  that  the  statistical  inversion  method  has  great 
potential  in  the  remote  determination  of  two  atmospheric  parameters,  the  struc¬ 
ture  constant  and  the  average  transverse  wind  velocity. 
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